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ABSTRACT 



Aims. After proposing a new method of deriving the atmospheric time constant from the speed of focus variations (Kellerer & Tokovinin 
2007), we now implement it with the new instrument, FADE. 

Methods. FADE uses a 36-cm Celestron telescope that is modified to transform stellar point images into a ring by increasing the central 
obstruction and combining defocus with spherical aberration. Sequences of images recorded with a fast CCD detector are processed to 
determine the defocus and its variations in time from the ring radii. The temporal structure function of the defocus is fitted with a model 
to derive the atmospheric seeing and time constant. We investigated by numerical simulation the data reduction algorithm and instrumental 
biases. Bias caused by instrumental effects, such as optical aberrations, detector noise, acquisition frequency, etc., is quantified. The ring image 
must be well-focused, i.e. must have a sufficiently sharp radial profile, otherwise, scintillation seriously affects the results. An acquisition 
frequency of 700 Hz appears adequate. 

Results. FADE was operated for 5 nights at the Cerro Tololo observatory in parallel with the regular site monitor. Reasonable agreement 
between the results from the two instruments has been obtained. 
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1. Introduction 

The site- and time-dependent performance of telescopes, and 
especially of interferometers, can be characterized by the pa- 
rameters seeing, so (or, equivalently, the Fried parameter ro = 
0.98/l/eo), and the coherence time, to, that determines the re- 
quired reaction speed of adaptive-optics (Roddier 1981). The 
variability of these parameters makes monitoring instruments 
essential. Seeing is usually measured with the Differential 
Image Motion Monitor, DIMM (Sarazin & Roddier 1 19901 1. 
However, a practical method of measuring to is still lacking. 
At present this parameter is inferred from the vertical pro- 
files of wind speed and turbulence, from the temporal anal- 
ysis of image motion, from scintillation, etc. (cf. the review 
in Kellerer & Tokovinin 2007 hereafter KT07). In particular, 
a Multi-Aperture Scintillation Sensor, MASS (Kornilov et al. 
2003 ) deduces the coherence time from scintillation, but this 
method (Tokovinin 2002) is only approximate and has not yet 
been verified by comparison with other techniques. 



A new method for measuring the coherence time with a 
small telescope has recently been proposed in KT07. This 
method, termed FADE (FAst DEfocus), is based on recording 
and processing focus fluctuations produced by the atmospheric 
turbulence in a small telescope. The variance of defocus is pro- 
portional to (D/ro) 5 ^ (Noll 1 1976b and gives a measure of the 
seeing, where D is the telescope diameter. As shown in KT07, 
the variance of the speed of defocus is related to to (this re- 
lation is given below in Sect. 13.41 where the derivation of to 
from a sequence of fast defocus measurement is explained). In 
principle, To can be obtained from the temporal analysis of al- 
most any quantity affected by turbulence. However, tilts, the 
easiest to measure, are typically corrupted by telescope shake 
and guiding errors, hence are not suitable. The DIMM instru- 
ment is immune to the wind shake, but it is intrinsically asym- 
metric. An early attempt to extract to from the DIMM signal 
by Lopez d 1 992b revealed the complexity of this approach and 
did not result in a practical instrument. If we discard tilts, the 
next second largest and slowest atmospheric terms are defocus 
and astigmatism. Defocus has angular symmetry and the rate 
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Fig. 1. Overview of the FADE instrument and data analysis. 



of its variation is related to to. Attractive features of the FADE 
method thus are: 

- direct measurement of To and bq, 

- use of the whole telescope aperture, 

- immunity to tilts, 

- small telescope size. 

FADE can be useful for site testing and monitoring, but its fea- 
sibility has so far been demonstrated only by numerical sim- 
ulation. Here we present an instrument implementing the new 
method. 

The instrumental set-up is described in Sect.[2] Section[3] 
outlines the data analysis algorithm, while various instrumen- 
tal effects are evaluated in Sect. [4] by numerical simulation. In 
Sect. [5] the seeing and coherence time measured with FADE 
are checked for consistency and are compared to simultaneous 
data from the DIMM and MASS instruments. Section|6] con- 
tains conclusions and an outline of further work. Mathematical 
derivations are given in AppendiceslAlandlBl 

2. The instrument 

2.1. Operational principle 

The temporal structure function (SF) of atmospherically- 
induced focus variations D/t(f) is related to the average wind 
speed in the atmosphere V. The initial, quadratic part of SF is 

D 4 (0 ~ Ct 2 f dh C 2 „(h) V 2 (h) = Ct 2 v\ f dh C 2 n (h), (1) 
Jo Jo 

where the proportionality coefficient C depends on the wave- 
length, telescope diameter and central obstruction (see KT07 
and Appendix B for the derivation), C 2 {h) and V(h) are verti- 
cal profiles of the refractive-index structure constant and wind 
speed. By measuring SF, we can estimate the integral in ([TJ that 
is similar to the integral entering the definition of to (Roddier 
[19811 

118^ 2 J dh C 2 „(h) V 5/3 (h)\ . (2) 

The exponent of the wind speed is slightly different; however, 
we show below that fitting the measured SF to a theoretical 
model leads to a good estimation of tq. 



The defocus aberration can be measured with a wave-front 
sensor of any type or can be simply inferred from the size of 
a slightly defocused long-exposure stellar image (Tokovinin 
& Heathcote 2006). For FADE, a simple, fast, and accurate 
method is required. We chose to introduce a conic aberration 
into the beam in order to form a ring-like image. A small de- 
focus slightly changes the radius of the ring. Ring-like images, 
"donuts", are obtained by defocusing a telescope with a central 
obstruction. However, unlike a donut, the ring is fairly sharp in 
the radial direction, which means that the determination of the 
ring radius is largely insensitive to intensity fluctuations (scin- 
tillation) at the telescope pupil. 

There is an inherent similarity between FADE and DIMM. 
In a DIMM, two peripheral beams are selected and are deviated 
by prisms to form an image of two spots. In FADE, the prisms 
are replaced by a conic aberration and the whole annular aper- 
ture is used to form a ring-like image instead of two spots. 

The ring images are recorded by a fast CCD detector and 
stored on a computer disk (Fig. [TJ. They are processed offline 
to determine a temporal sequence of ring radii, pit), related to 
the defocus. In order to estimate the atmospheric parameters so 
and To, the temporal structure function of the defocus variations 
is computed (Eq. [T4T i and fitted to a model. 

Atmospheric defocus fluctuations are fast: their temporal 
correlations decrease with a half-width 0.3 times the aperture 
crossing time f closs = D/V, i.e. with 2.2ms for a telescope di- 
ameter D = 0.36 m and wind speed V = 50m/s (cf. Appendix 
B). To capture the focus variations of interest, an acquisition 
frequency v > 500 Hz (exposure time < 2 ms) is required, 
which is attainable with today's fast CCD detectors. 

2.2. Hardware 

We assembled the FADE prototype from readily available com- 
mercial components (Tabl^TJi. A 36-cm telescope was selected. 
In a smaller telescope, the focus variations are smaller, (a 2 ) cc 
D 5 ^ 3 , and faster, f closs cc D, hence more difficult to measure. 
Use of a fast CCD - GE 680 from Prosilica - is critical for the 
instrument, because it permits continuous acquisition with an 
image frequency 740 Hz when a 100x100 subsection of the full 
frame (region-of-interest, ROI) is read out. The signal is digi- 
tized in 12 bits. With the lowest internal gain setting, dB, the 
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Table 1. Components of the FADE instrument 



Component 



Description 



Telescope 
Central obstruction 
Aberrator 

Detector 

Interface 
Computer & OS 



Celestron C14, D = 0.356 m, F = 3.910 m 
Circular mask of 150 mm diameter 
PCX lenses (Linos 312321 & 314321), 
d L = 25 mm, f L = ±50 mm 
Prosilica GE680, 640 x 480, 
pixel 7.4 fim (0.39") 
Gigabit Ethernet IEEE 802.3 lOOObaseT 
De/ZD410, Windows XP 



conversion factor 2.86 ADU per electron and the readout noise 
38ADU = 13.4e were measured. According to the specifica- 
tions, the maximum quantum efficiency (QE) is 0.5 electrons 
per photon at wavelength A = 0.50 fim with a full-width half- 
maximum spectral response of roughly AA = 0.25 yum. Indeed, 
the measured fluxes from stars correspond to the overall system 
QE of 0.35-0.40, including atmospheric and optical losses. 



2.3. Optics 




Fig. 2. Optical layout of FADE (not to scale). The stellar light 
comes from left to right. After reflection from the telescope pri- 
mary and secondary mirrors Ml and M2, the converging beam 
passes through the "aberrator" consisting of the lenses LI and 
L2 and is detected by the CCD. The wave-front cross-sections 
before and after the aberrator are illustrated by thick line seg- 
ments. 



To create annular images, a conic aberration must be intro- 
duced into the beam. Conic lenses, axicons, have wide tech- 
nical and research applications and are commercially avail- 
able. In the case of FADE, we can successfully approximate 
a conic wave-front by a combination of a quadratic (defocus) 
and higher-order (spherical aberration) terms. 

An aperture with a relative central obstruction e gives a 
ring-like defocused image with the average angular radius p = 
V48/[D(1 - e)], where «4 is the defocus coefficient in lin- 
ear units (throughout this article, the Zernike aberration coeffi- 
cients are given in the Noll 1 1 9761 notation) . Such a ring can be 
sharpened by adding spherical aberration. Elementary analysis 
shows that for a relative central obstruction e ~ 0.4, the devi- 
ation from a conic surface is minimized when an = -0.1 a\. 
The rms error of this approximation is 0.01 04. In order to get a 



ring radius of 5" with a 35-cm telescope, we need the defocus 
amplitude 04 m 600 nm; hence, the difference of the approxi- 
mated wave-front from a perfect cone will be 6 nm, or A/80 at 
A = 500 nm. 

To attain the desired aberration, we used an assembly of 
two simple piano-spherical lenses with equal but opposite cur- 
vature radii, which can be seen as a plane-parallel plate con- 
taining a meniscus-shaped void (Fig. |2). The thickness of the 
meniscus is adjusted by changing the gap g between the lenses. 
The positive lens is closer to the primary mirror, so that the 
meniscus curvature opposes the curvature of the wavefront. We 
used lenses with focal lengths fa = ±50 mm and a gap g = 
0.7 mm. When this element is placed at distance / = 93.5 mm 
in front of the detector and the telescope is suitably refocused, 
a ring image of radius p * 5" is formed. Optical modeling 
in Zemax shows that this "aberrator" is reasonably achromatic. 
The spherical aberration is proportional to g I 6 , therefore it can 
be adjusted over a wide range. To block the inner part of the 
wavefront where it deviates from the cone, a central obstruc- 
tion of 150 mm diameter, i.e. a relative diameter e = 0.42, 
was placed at the telescope entrance. Obviously, there are many 
other possible optical arrangements to obtain wave-fronts with 
spherical aberration. 

The average ring image in the real FADE instrument 
(Fig. [3) shows marked aberrations other than conical, caused 
by the defects of optical surfaces and of alignment. Similar 
rings were reproduced in our simulations with a combination 
of coma and higher-order aberrations (cf. Sect. 14.2b . We also 
fitted the Zernike aberrations directly using the donut method 
(Tokovinin & Heathcote 2006) and found that the coma coef- 
ficient could reach ~100nm (1.2rad). Furthermore, the defo- 
cus was not always kept at its optimum value as required for 
sharp ring images. The effect of such aberrations is studied in 
Sect. l4.4l bv simulation. 

2.4. Acquisition software 

Since the GE680 detector is relatively new, with no readily 
available software development kits as yet, we used the com- 
mercial software, Streampix from Prosilica. It provides all nec- 
essary functions for detector control and data storage in the 
FITS format, but the parameters need to be set manually at each 
acquisition, which requires constant attention. And they are not 
logged into the FITS headers or otherwise. Thus, Streampix is 
only a temporary solution. We checked that the image sequence 
is acquired at regular intervals, without time jitter. The detector 
was exposed for this purpose to a strictly periodic light signal 
at 10 Hz, and a 100 x 100 ROI was read at 400 Hz. The power 
spectrum of the flux calculated from these data is a narrow peak 
at (10.0 ± 0.2) Hz without significant tails. 

2.5. Observations 

The FADE instrument was installed in the USNO dome of 
the Cerro Tololo Inter- American Observatory (CTIO) in Chile 
for the period October 27 to November 3, 2006. The instru- 
ment was elevated ~4m above the ground. We pointed FADE 
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Fig. 3. From left to right: Simulated ring image - Image of Sirius - Average of 1024 simulated images - Average of 1024 Sirius 
images. The sequence of Sirius images was recorded on Nov. 2 at 6:46 UT. The parameters for the data and simulations are given 
in Sect.l4T2l 



at bright stars, Fomalhaut (a PsA, A3V, my = 1.16) in the 
evening, then Sirius (a CMa, A1V, my = -1.47). The expo- 
sure time ranged from 1 ms to 1 .9 ms for Fomalhaut and was 
< 0.5 ms for Sirius to avoid saturation. Figure [3] shows typical 
instantaneous and average images of Sirius, as well as simu- 
lated images. During our test run, the seeing was not very good 
at roughly 1", and the turbulence in the high atmosphere was 
strong and fast, as shown by the MASS data. 

3. Data analysis 

A correct algorithm of data processing and interpretation is 
critical for deriving the atmospheric parameters sq and to. We 
carefully selected the most robust method of calculating atmo- 
spheric defocus from the ring-like images and used numerical 
simulation to study the influence on the results of various in- 
strumental effects and of optical propagation (Sect.|4]below). 

3.1. Estimating the ring radius 

The center of gravity of the image (x c ,y c ) is calculated by the 
usual formula 

Xc = 2j Xl ' k 2j Il ' k and y ° = 2j y, ' k Il ' k l 2j /u " ^ 

l,k l,k l,k l,k 

The ring radius p can then be estimated in a similar way, as the 
intensity-averaged distance from this center: 

p=Yj r '* i ^Yj Ii - k - (4) 

l,k l,k 

Here is the light intensity at pixel (/, k), and is the dis- 
tance of this pixel from the center. There are various caveats 
below the apparent simplicity of this procedure. 

There is no unambiguous way to assign a center to a real 
(distorted and noisy) ring image. A simple center-of-gravity is 
a very rough estimate of (x c ,y c ); in particular, it is affected by 
the intensity fluctuations in the ring due to scintillation. It is 
better to compute (x c ,y c ) with clipped intensities: below a 
threshold and 1 above, the threshold being set safely above the 
background and its fluctuations. This initial estimate can be im- 
proved further by minimizing the intensity-weighted mean dis- 
tance of the pixels from the ring, as described in Appendix A. 
However, small inaccuracies in the center determination do not 



affect the resulting radius critically, and in fact, we found the 
initial estimate to be adequate. 

A second caveat concerns the choice of the pixels used for 
the radius estimate. A considerable fraction of pixels lie outside 
the ring in an empty area that only contributes noise. To reduce 
the noise with a minimal loss of information, we restricted the 
pixels used in © to a mask of inner radius p - Ap and outer 
radius p+ Ap, where p is the average ring radius. We express the 
mask half-width Ap as a fraction 6 of the diffraction half-width 
of the ring, 

Ap = 6 /l/[0.5 D (1 - e)]. (5) 

Figured shows that a mask with 6-2 would be good for 
an ideal, diffraction-limited ring. For a typical image sequence, 
however, the ring is widened by telescope aberrations and at- 
mospheric distortions, so we set 5 = 4, which covers the actual 
ring image with a sufficiently conservative, but still reasonable, 
margin. 

The simulations show that scintillation and aberrations add 
to the fluctuations of the estimated radii and thus bias the results 
of FADE. To reduce this effect, we subdivide the ring into eight 
45° sectors and - utilizing the same center estimate (x c ,y c ) - 
apply Eq.|4]to each sector separately, and then average the re- 
sult. This reduces the effect of azimuthal intensity variations. 
An added advantage of the procedure is that the relative vari- 
ance s of the total intensities in the sectors 4 with respect to 
their average Ik serves as a measure of the scintillation, hence 
of the turbulence height, 

1 8 

s = - Yjih ~ h) 2 I h 2 . (6) 

The method of calculating the ring parameters (x c ,y c ,p) is 
less rigorous than fitting a wave-front model directly to the im- 
age. The big advantage of the estimator however, is its sim- 
plicity. 

3.2. Noise and limiting stellar magnitude 

The errors of the radius estimates caused by photon and read- 
out noise are obtained by differentiating Eq. © and using the 
independence of the noise in each pixel: 
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Fig. 4. Total intensity inside concentric circles of radii p for 
the average of 1024 centered images. Full line: sequence of 
Fomalhaut images recorded on Nov. 2. Dashed line: simulated 
diffraction-limited ring images (see Table|2]). 



Recall that the atmsopheric defocus a$ is related to the phase 
distortion ip(r) as 



fl4 



-/ 



d 2 r z 4 (r) ip(r), 



(9) 



where r is the normalized coordinate vector on the pupil and 
Z4(r) the orto-normal Zernike defocus given by Noll ( 1976 ) for 
the circular aperture and in Fig. |5]for the annular aperture. 
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Fig. 5. Response functions z(r) on annular aperture for Zernike 
defocus, conic aberration, and average radial gradient. The first 
two functions are normalized in the Noll (119761 1 sense. The co- 
efficients «4, a e , and a g are calculated as integrals ([9). Here 6 is 
the Dirac's delta function. 



S 2 P = Yjhk(r,, k -pf 

l,k 



ph- 



(8) 



Here, cr ron is the rms detector noise, A^ p h = 2/,* h,k the total stel- 
lar flux in one exposure (both in electrons), p the average ring 
radius, and is the distance of pixel (/, k) from the center, ex- 
pressed either in pixels or arc-seconds. The rms ring-width 5 P 
quantifies the ring sharpness, which turns out to be critical for 
getting unbiased measurements with FADE (see Sect. 14.4b . The 
summation is extended only over pixels inside the mask, as de- 
scribed in Sect. 13.11 We recognize a familiar sum of the readout 
noise (first term) and photon noise (second term), where the 
first term typically dominates. Equation (0 does not account 
for such additional noise sources as scintillation, image distor- 
tion, etc. 

Formula (0 is useful for predicting the limiting magnitude 
of FADE. A star of zero V-magnitude gives a flux A^ p h ~ 6 • 10 5 
photo-electrons in 1 ms exposure in our instrument. The rms 
noise on the radius estimate with plausible parameters (p = 5", 
0"ron = 13, 6 = 4) is then about 2 mas. It will increase to 
20 mas for a star with my = 2.5"' - still much less than the 
atmospheric signal. Despite very short exposures, FADE is not 
photon-starved. 

3.3. The response coefficient of FADE 

The relation between the ring radius fluctuations Ap and the 
atmospheric defocus (Zernike coefficient as) is intuitively clear. 
But what is the exact coefficient A in the formula Ap = Aa\l 



A reaction of our simple radius estimator © to a small per- 
turbation of phase and amplitude at the telescope pupil can be 
determined analytically (cf. Perrin et al. l2003l for an example of 
similar analytics). It turns out that the response to a phase per- 
turbation in the pupil plane is not exactly proportional to the 
Zernike defocus. Moreover, it depends on the adopted mask 
half- width 6. For 5 ~ 1, the response resembles a cone, so that 
FADE measures something similar to a conic aberration. On 
the other hand, for 6 > 2 the computed ring radius is related 
to the average radial gradient of the wave-front, and therefore 
FADE measures the difference a g between the phase averaged 
on the outer and inner edges of its annular aperture. Its response 
is further modified when the ring is distorted by aberrations. 
In this case, the radius estimate is sensitive to both amplitude 
and phase fluctuations. Although we developed a full analyti- 
cal treatment of this problem, it is omitted here for the sake of 
simplicity. 

The three quantities - Zernike defocus a\, conic aberration 
a c , and average phase gradient a g - are similar, especially on 
the annular aperture (Fig. [Sj. FADE measures something else, 
but its response is most closely approximated by a g when the 
ring radius is calculated with a large mask width 5. Let a g be the 
average phase difference between the outer and inner borders 
of the aperture, the corresponding change of the angular ring 
radius is then 



nD{\ - e) 



(10) 



The Zernike defocus on the annular aperture is proportional to 
04 VT2r 2 /(l - e 2 ), where r is normalized by the pupil radius. 
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Hence, a g — a 4 x Vl2, and the proportionality coefficient A 
follows from Eq. (TTOb , 



Ap = A a 4 — a 4 — 



VT2 



ttD 1 -e 



(11) 



The atmospheric variance of the defocus a 4 or gradient a g on an 
annular aperture can be computed, as done by Noll (119761 1 for 
a filled aperture. Alternatively, the variance of the ring radius 
may be directly written as 



al = C p (A/D) 2 (D/r Q ) 



5/3 



(12) 



analogous to similar formulae for the gradient or Zernike tilt. 
Our numerical calculation for the average-gradient response 
(see Eq.fTOt leads to an approximation valid for e < 0.6 with 
an accuracy of +7 10~ 5 : 



C p ~ 0.03288 + 0.0503e - 0.05638e 2 + 0.04056e 3 . 



(13) 



We studied the response of FADE by analytical calculation 
and numerical simulations and found that the exact coefficient 
C p in Eq. [12] depends on all parameters of the instrument and 
data processing. A choice of 6 > 2.5 ensures a relative stability 
of the response with respect to small aberration, propagation, 
etc. A small correction to the "ideal" response coefficient is 
finally determined by simulation (Sect. 14. 3t and applied to the 
real data. 

The lack of a unique, well-established coefficient relating 
measurements to atmospheric parameters may appear disturb- 
ing. However, a similar analysis applied to the classical DIMM 
instrument leads to the conclusion that its response, too, de- 
pends on the details of centroid calculation and, furthermore, 
is modified by propagation and optical aberrations. In this re- 
spect, FADE and DIMM are not different. 

3.4. Derivation of the seeing and coherence time 

We convert the measured ring radius into defocus using coeffi- 
cient A (see Eq. fTTb and calculate the temporal structure func- 
tion of defocus D 4 (t), 



D 4 (t) = <[a 4 (f' + f) - a 4 (f')] 2 >- 



(14) 



A typical SF is plotted in Fig. [6] A theoretical expression for the 
defocus SF has been derived in KT07. We generalize it to annu- 
lar apertures in Appendix B. The initial, quadratic part of SF is 
directly related to the combined time constant tq of all turbulent 
layers. However, the acquisition frequency is not fast enough to 
capture the initial quadratic part of SF extending only to time 
lags of < 0.1f cross . In order to get two points on this part for a 
layer moving with V = 50 m/s, a frame rate of ~ 3 kHz would 
be required. 

To overcome the sampling problem, we fit the initial part of 
SF to a model of N turbulent layers with Fried parameters roj 
and velocities V,, 1 < i < N: 



D 4 (t > 0) = 1.94D 5/3 r Q ^ 3 K 4 (2tVi/D, e) + 



i=l 



2cr 2 . 

p, noise 

A 2 ' 



(15) 



where the function K 4 (J3, e) is defined in Appendix B and 
cr 2 is the noise of the radius estimate determined by Eq. [7] 

p, noise J " — 1 
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Fig. 6. Structure function of focus variations measured at 
700 Hz (crosses) fitted with a model of three turbulent layers 
(line). 



The adjusted parameters are ro,,- and V,. As will be seen in 
Sect. |5.1| the estimate of to is independent of N if N > 3. 
Accordingly, a three-layer model is chosen for the data anal- 
ysis (Fig. |6j. We fit only the initial part of SF, up to the 
time increment At. Its exact value is not critical, as long as 
it is large enough for unambiguous fitting of the parameters, 
At v > 2N + 1 . For further data analysis, we set At = 40 ms. 
The atmospheric parameters (ro, V, to) are calculated as 



r Zj r 0,i 



-5/3 



i=l 
N 

(V/rof 3 = Yj( V t/ r 0,if 3 < 

!=1 

T = 0.314 r /V. 



(16) 

(17) 
(18) 



The estimate of ro is also obtained directly from the ring-radius 
variance cr 2 by subtracting the noise, 



p,noise 



CJA/D)\D/r ) 



5/3 



(19) 



When SF reaches its asymptotic value on time increments 
smaller than 40 ms, the same value of ro is derived from the 
ring-radius variance (Eq.[T9ll and from the model (Eq.fTBTl. The 
robustness of parameter estimates derived by model fitting has 
been confirmed by numerical simulation (Sect.|4]i and by fitting 
alternative models to real data (Sect. 15. It . 

4. Simulations 

A new seeing monitor can be validated by comparing it with an- 
other, well-established instrument. In the case of FADE, how- 
ever, there is no reliable comparison data on tq. Instead, we 
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Table 2. Simulation parameters. 
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simulated our instrument numerically as faithfully as we could 
and studied the influence of various instrumental and data- 
reduction parameters on the final result. 

4.1. Simulation tool 

Our simulation tool generates the complex amplitude of the 
light field propagated through one or several phase screens with 
Kolmogorov spectrum. The screens are typically 1024 2 pixels 
with 1 cm sampling, i.e. about 10 m across. The resulting am- 
plitude pattern is periodic, without edge effects. It is "dragged" 
in front of the simulated telescope with a chosen wind speed, 
wrapping around edges in both coordinates and eventually cov- 
ering the whole area. The monochromatic images created by a 
telescope with a perfect conic aberration of specified ampli- 
tude and, possibly, some additional intrinsic aberrations are re- 
binned into the detector pixels, distorted by readout and pho- 
ton noise and fed to the data-analysis routine instead of the real 
data. Our tool has been verified by comparing with analytical 
results for weak perturbations and has been used for simulating 
other instruments such as DIMM and MASS. The limitations 
of this tool are the monochromatic light, single wind velocity 
for all layers, and instantaneous exposure time. 

We used two alternative, nearly equivalent ways of produc- 
ing ring images. In the first method, a perfect conic wavefront 
was generated, and its amplitude was expressed as a ring ra- 
dius p. In the second method, we do not apply conic aberra- 
tion (p = 0), but instead select a combination of defocus and 
spherical aberrations to mimic a real telescope. The sense of 
the Zernike coefficients a\ and an in both cases is distinct. 

4.2. Parameters of the simulations 

For convenience, the simulation parameters are gathered in 
Table 12 These parameters are: 

- number of images in the sequence N, 

- acquisition frequency v, 

- exposure time df, 

- visual stellar magnitude m, 

- readout noise cr ron , 

- conic aberration quantified by the average ring radius p, 

- amplitudes of the Zernike aberrations a\ (defocus), a-i 
(coma), an (spherical), and 027, 



- altitude of the single turbulent layer h, 

- seeing sq, 

- wind speed V. 

Simulated ring images are compared in Fig.[3]to the images 
of Sirius recorded on Nov. 2. The combination of exposure time 
and magnitude results in the detected flux of 3 • 10 5 electrons 
per simulated image, as in the actual images of Sirius. For 
this sequence, the estimated turbulence parameters equal: sq = 
1.05", V = 17 m/s, To = 1.86 ms (Fig. [6]). The same parameters 
are chosen for the simulated images. For the closest resemb- 
lence between simulated and real images, telescope aberrations 
are set to a-j = 0.7 rad, a\i = -0.75 rad, an = 0.3 rad. The tur- 
bulence is placed at 13 km altitude to reproduce the actual level 
of scintillation, evaluated from the intensity variance between 
ring sectors s = 0.011 (cf. Eq.|6]l. 

4.3. Refining the response coefficient 

The coefficient A relating radius variation to defocus is given 
by Eq. [TT] Its actual numerical value, however, depends on the 
method of radius estimation and, in particular, on the choice 
of the mask width 6. For 6 — 4, we determined it to equal 
1 .077 by comparing to estimates from sequences of simulated 
images, to the nominal input value of To when Eq. Q~T]is use d 
to relate the radius to defocus. The corresponding simulation 
parameters are summarized in Table[2] Thus, 

Ap/a 4 = 1.077 A 7"— ■ (20) 
nD 1—e 

4.4. Instrumental biases 

The data analysis relies on radius estimates that can be al- 
tered by telescope aberrations, scintillation, detector and pho- 
ton noise, etc. Here we evaluate the instrumental bias by chang- 
ing some parameters, while other parameters are fixed. In each 
case, a sequence of 1024 simulated images is generated with 
the parameter values listed in Table|2] The wind speed is set to 
35 m/s and the coherence time is then changed by modifying 
the seeing. 

Ring sharpness. The ring is sharp in the radial direction when 
the wavefront is exactly conic. A good approximation of the 
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Fig. 7. Influence of the ring sharpness - quantified in terms of the ring width S p - on the seeing and coherence time estimates. 
Top - turbulence layer at ground level, bottom - turbulence layer at 5 km altitude. Simulation parameters are given in Table|2] 



conic wavefront is achieved by the optimum combination of 
defocus and spherical aberrations, a u = -0.1 a^. Here we ex- 
plore the effect of unsharp ring images by setting 04 = 12rad 
and varying a 1 1 about its optimum value an = - 1. 2 rad. Unlike 
the rest of the simulations, we do not apply conic aberration and 
set p = 0. Wrong values of an make the ring wider, as shown 
by its rms width, 5 P (Eq.0. 

As seen in Fig. [7] the seeing and coherence-time estimates 
are biased in the case of blurred rings and high-altitude turbu- 
lence. The sign of the bias depends on the sign of the deviation 
from the optimum ci\\. When the turbulent layers are low, the 
scintillation is weak and the parameters are correctly derived 
even if the ring images are blurred. 

To ensure a correct derivation under any atmospheric con- 
ditions, the ring width should be close to its diffraction-limited 
value 6 p fi'. 

6 P < l.2S pJ0 ; 8 p ,o = 1.7A/[D(l-e)], (21) 



where the coefficient 1.7 is determined from the width, 6 P = 
0.87", of diffraction-limited rings. Given the instrumental set- 
up, images should be rejected if 5 P < I". However, all images 
recorded with the FADE prototype have 5 P > 1". Hence, we 
apply a softer data-selection criterion: 5 P < 1 .25", and note that 
the resulting estimates might still be biased if the turbulence 
was high. 



Stellar magnitude, ring radius, and coma. Figure[8]examines 
the stability of the seeing and coherence time estimates with 
respect to the stellar magnitude my, ring image radius p, and 
coma aberration a-/. In agreement with Eq.|7l estimates are cor- 
rect up to stellar magnitudes 2-3. Spatial sampling and coma 
aberration do not affect the estimates if p > 2" (i.e. 5 pixels) 
and a-/ < 2 rad. 
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Fig. 8. Dependence of the seeing and coherence time estimates on stellar magnitude my, mean ring radius p, and coma aberration 
a-j. Simulation parameters are given in Table [2] 



5. Analysis of observations 

Seeing and coherence time were estimated from all sequences 
of 4000 images recorded with FADE at Cerro Tololo between 
October 29 and November 2, 2006. In this section, we check 
the FADE results for consistency and compare them with the 
MASS-DIMM. 

5.1. Influence of instrumental parameters 

During data acquisition, instrumental parameters were varied 
over a broad range to evaluate their effect on the results. Even 
though the non-stationarity of the atmosphere precludes direct 
comparisons, some conclusions can nevertheless be drawn. 

The ring sharpness, 5 P , has been identified as a major 
source of instrumental bias in FADE when significant high- 
altitude turbulence is present. In our data, most images have 
1" < dp < 1.5", whereas a perfect diffraction-limited ring 
has 6 P = 0.9". Analysis of the average ring images confirms 
that the optimum combination of defocus and spherical aberra- 
tions was not reached, a\ and an often having the same sign 
rather than opposite signs. For our data, the dispersion and 
mean of to increase when 6 > 1.25"; accordingly, sequences 
with 6 P > 1.25" are disregarded. Still, some bias caused by 
radially defocused images remains. 

As described in Sect. 13. 41 the data are fitted to a model with 
a discrete number of turbulent layers, N. What is the minimum 
value of N that permits a correct derivation of to? Figure|9] 



t 

c 

t- 




0.5 1.0 1.5 2.0 2.5 3.0 
t , ms - 3 layers 



Fig. 9. Coherence time derived from the recorded data by fit- 
ting a model with N - 3 layers (x-axis) compared to the coher- 
ence time derived with N - 2 (dotted line) and N - A (solid 
line) models. 



shows that the To values obtained with 3 and 2 (resp. 4) lay- 
ers differ on average by 7% (resp. 2%). An average difference 
of 2% is likewise obtained when comparing the estimates with 
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3 and 5 layers. A 3-layer model is thus a good compromise 
enabling a fit to the data with only six parameters. 




Fig. 10. Coherence time derived from all sequences recorded 
with v > 700 Hz when every image (x axis) and every other 
image (y axis) is considered. 



The influence of the acquisition frequency on the measured 
coherence time is examined in Fig.[l0] The data sequences 
recorded at frequencies v > 700 Hz were re-analyzed consid- 
ering every other image. The coherence time obtained with a 
slower v/2 sampling is on average 9% longer than with the 
fast sampling. This difference is reproduced by simulations if 
the turbulence is placed at 5 km altitude and if the ring-images 
are slightly defocused in the radial direction (a\\la$ ~ -0.07 
or -0.14 instead of 011/04 = -0.11 corresponding to a sharp 
ring). The effect of temporal under-sampling is perceptible if 
the same comparison is repeated with sequences recorded at 
frequencies below 700 Hz: the number of points on the initial, 
increasing part of SF is then not always sufficient to unambigu- 
ously extract the six fitted parameters, and the coherence time is 
poorly constrained. To ensure a correct temporal sampling un- 
der fast turbulence, we ignored the sequences with v < 500 Hz. 

In line with the simulations, the coherence time estimates 
do not depend on the average ring image radius. Similarly, the 
parameter statistics seem unbiased by the stellar flux and by the 
exposure time. While the sequences of Sirius (my = —1.5) and 
Fomalhaut (my = 1.2) images were recorded with exposure 
times of dt < 0.5 ms and 1.0 < df < 1.9 ms, respectively, no 
obvious difference exists between the mean and rms of the 
atmospheric parameters measured in terms of these two stars, 
4 ■ (0.9 + 0.2)" el : (0.8 ±0.1)" 

Tq : (1.4 + 0.5) ms 7*: (1.3 ±0.5) ms. 



5.2. Comparison with MASS and DIMM 

In this section, the seeing and coherence time obtained with 
FADE are compared to simultaneous measurements by the 
CTIO site monitor located at 10 m distance from FADE on a 
6 m high tower. The monitor consists of a combined MASS- 
DIMM instrument fed by the 25-cm Meade telescope and look- 
ing at bright (V = TP 1 .. .3™) stars near zenith. Of particular in- 
terest here is the time constant t estimated by MASS from 
the temporal characteristics of scintillation by the method of 
Tokovinin (2002). This method is intrinsically biased because it 
does not account for the turbulence below ~ 500 m. Moreover, 
it has been recently established by simulations that the coeffi- 
cient used to calculate to in the MASS software must be in- 
cr eased by 1.270 In the following, we correct to by applying 
this coefficient and including the contribution of the ground 
layer: 

t 5/3 = (1.27 T M ASsr 5/3 + 1 18 A' 2 V 5 ^ (Cl d/i) GL . (22) 

The turbulence integral in the ground layer, (C 2 d/i)cL, is com- 
puted from the difference between the turbulence integrals 
measured by DIMM (whole atmosphere) and MASS (above 
500 m), while the ground layer wind speed, Vql, is known 
from the local meteorological station. Even after correction 
by Eq. (T22l . the coherence time measured by MASS-DIMM 
should be taken with some reservation because it has never 
been checked against independent instruments and some bias 
is possible. 

FigureQj] compares the estimates of eq and to obtained 
with FADE from October 29 to November 2 to the results of 
MASS-DIMM. We do not expect detailed correlation, because 
the instruments were sampling different atmospheric volumes. 
As seen in Fig.Qj] the seeing measurements are better corre- 
lated than the coherence times. 

Statistically, it appears that FADE slightly underestimates 
the seeing. This effect is reproduced with simulations of high- 
altitude turbulence if the ratio of spherical aberration to defocus 
On /a4 is set higher than its optimum value -0. 1 corresponding 
to sharp ring images. In this case, FADE also underestimates 
the coherence time. The bias on to can, however, not be as- 
certained by Fig.QT]because the T o estimates by MASS might 
likewise be biased. Thus, the comparison presented in this sec- 
tion cannot be considered as a validation of FADE. 

6. Conclusions and perspectives 

We have built a first prototype of the site-testing monitor, 
FADE, suitable for routine measurements of the atmospheric 
coherence time to, as well as the seeing eo. The instrument 
was tested on the sky. Extensive simulations substantiate the 
validity of the FADE results and indicate potential instrumen- 
tal biases. Our main conclusions are as follows: 

- The sampling time of the image sequence must be a small 
fraction of the aperture crossing time f cross = D/V (~10 ms 

1 See the unpublished report by Tokovinin (2006) at 
http://www.ctio.noao.edu/~atokovin/profiler/timeconst.pdf 
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Fig. 11. Seeing and coherence time measured with FADE be- 
tween October 29 and November 2, compared to simultaneous 
measurements by the MASS-DIMM. The average values and 
standard deviations of parameters and the correlation coeffi- 
cients are indicated. 



for D = 0.36 m and wind speed V = 36 m/s). Sampling at 
v > 500 Hz appears adequate under most conditions. 

- The sharpness of the ring image in the radial direction does 
not bias the results when the turbulence is located near the 
ground. But it can bias both to and so estimates when high 
layers dominate, and strict control of the telescope aberra- 
tions is thus required. The aberrations (hence the data va- 
lidity) can be evaluated a posteriori from the average ring 
image. Real-time estimates of the ring radius p and width 
5 P are needed to ensure good optical adjustment of the in- 
strument. 

- The FADE monitor with 36-cm telescope can work on stars 
as faint as m v = 3 m . 



- A simple estimator of the ring radius (Eq.|U) is adequate 
and robust, provided a wide enough mask around the ring 
{5 ~ 4) is used in the calculation. 

- Moderate telescope aberrations such as coma are accept- 
able. The results are not critically influenced by small tele- 
scope focus errors. 

The current FADE prototype stores all image sequences, 
leading to a large data volume; the data are processed offline. 
While this procedure was necessary for the first experiments, 
online processing will be implemented in a definitive instru- 
ment. We have formulated and tested the data processing algo- 
rithm and can now develop adequate real-time software. 

We plan to develop an improved version of FADE with real- 
time data analysis. It will be compared to simultaneous esti- 
mates of the atmospheric time constant from currently work- 
ing adaptive-optics systems (Fusco et al. 2004) and/or long- 
baseline interferometers such as VLTI. Characterization of 
Antarctic sites for future interferometers is an obvious appli- 
cation for FADE. 
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Appendix A: Estimator of the ring radius and 
center 

The parameters of the ring-like image - its center (x c , y c ) and 
radius p - can be derived by minimizing the intensity- weighted 
mean squared distance of the pixels from the circle, 5 2 : 



= J]hk(n. k -P) 2 1 J]i,, k , 

It, Lk 



(A.1) 



where = [(/ - x c ) 2 + (k - y c ) 2 ] ' 5 is the distance of pixel 
(/, k) from the ring center, (x c ,y c ). Setting the partial derivative 
of S 2 over p to zero, we obtain the radius estimator of Eq.@] 
However, it still depends on the unknown parameters (x c ,y c ). 
By use of Eq.[4] Eq. IA.ll is simplified to: 



2 E« htfu 
6 p = 



Tji,k hk njt 



(A.2) 



This formula does not contain p. The center coordinates (x c , y c ) 
can be derived by setting the partial derivatives of 5 2 p over pa- 
rameters to zero and solving the equations. We determine the 
center numerically by minimizing Eq. lA.2l and using the center- 
of-gravity coordinates as a starting point. 



Appendix B: Structure function of atmospheric 
defocus 

The temporal structure function of atmospherically-induced 
defocus variations - Zernike coefficient a^ in Noll's ( 1976) no- 
tation - has been derived in KT07 for a filled circular aperture. 
Here we generalize it to an annular aperture. Without repeating 
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the whole derivation, we refer the reader to KT07 and mod- 
ify only the spatial spectrum of the Zernike defocus, taking the 
central obstruction ratio e into account. The resulting expres- 
sion is 

D 4 (t) = 0.821 k 2 D 5/3 Ah C n (h) 2 K 4 e| ,(B.l) 

K 4 Q3,e) = j-j dx x- s/3 [I - J Q Q3x)] 

(1 - e 2 ) 4 Jo 



J 3 (x) A h{ex) i ^/iW 



(B.2) 



where £ = 27r//i, /„ is the Bessel function of order n, C n {h) 2 and 
V(/i) are the altitude profiles of the refractive-index structure 
constant and wind speed, respectively. Considering the known 
relation between the turbulence integral and the Fried param- 
eter, r 5/3 = 0.423 k 2 C\dh, we can also write the defocus SF 
produced by a single layer as 

D 4 (t) = l.94(D/r ) 5/3 K 4 (2tV/D,e). (B.3) 

For calculating the function K 4 , it is convenient to approximate 
the integral (IB.21 > by an analytical formula, as in KT07. We 
suggest the approximation 



K 4 (j3,e) 



1 + C 3 /3° + /J 6 ' 
where the coefficients are cubic polynomials of e: 



(B.4) 



equivalently, the parameter C\ in Eq. IB.4l ) is reduced by 0.8 for 
an exposure time f exp ~ 0.3D/V and a layer moving with the 
speed V . To adequately sample the SF features produced by the 
fastest-moving layers, the sampling time (hence exposure time) 
must be shorter than Q.3D/ V max , so the bias caused by the finite 
exposure in FADE can be neglected. In hindsight, this result 
could be expected: to follow the focus variations, we need such 
a fast sampling that the integration during the sampling period 
has a negligible effect. 
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C; - C,,o V £ k , 



(B.5) 



k=o 



cf. Table IB. 11 This approximation is valid for e < 0.6 with 
a maximum relative error of less than 5% (3% for e = 0.42) 
and correct asymptotes. The asymptotic value K 4 (oo, e) = C2 
gives the focus variance on annular aperture, analogous to the 
Noll's coefficient. For e — 0, we get C2 = 0.024 and the focus 
variance coefficient of 1.94 x 0.024/2 = 0.0233, in agreement 
with Noll's result. 



Table B.l. Coefficients of dB4 



Param. 


Co 


6 U 


e 1 


2 

e~ 


6 j 


c, 


0.04642 


1 


-0.182 


-2.431 


2.028 




0.0240 


1 


-0.017 


-3.619 


2.833 


c 3 


1 


1.25 








7.5 


a 


1 


2.18 


-0.93 









The function K 4 (J3, e) reaches half its saturation value at 
P = 0.63; hence, the atmospheric defocus correlation time is 
~ 0.3D/V, as is well known in adaptive optics. 

The above analysis is valid for instantaneous measure- 
ments, while the defocus is in fact averaged over the exposure 
time. This effect is usually non-negligible for the DIMM. The 
time averaging can be included as an additional factor in the 
integral (IB.2l i. as done e.g. in Tokovinin (2002). We made this 
calculation and found that the initial, quadratic part of D 4 ( f) (or, 



